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We show that, for an exactly solvable quantum spin model, a discontinuity in the first derivative 
of the ground state concurrence appears in the absence of quantum phase transition. It is opposed 
to the popular belief that the non-analyticity property of entanglement (ground state concurrence) 
can be used to determine quantum phase transitions. We further point out that the analyticity 
property of the ground state concurrence in general can be more intricate than that of the ground 
state energy. Thus there is no one-to-one correspondence between quantum phase transitions and 
the non-analyticity property of the concurrence. Moreover, we show that the von Neumann entropy, 
as another measure of entanglement, can not reveal quantum phase transition in the present model. 
Therefore, in order to link with quantum phase transitions, some other measures of entanglement 
are needed. 
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Quantum entanglement, as one of the most fascinat- 
ing feature of quantum theory, has attracted much at- 
tention over the past decade, mostly because its non- 
local connotation Q is regarded as a valuable resource 
in quantum communication and information process- 
ing [2(. Recently a great deal of effort has been de- 
voted to the understanding of the connection between 
quantum entanglement and quantum phase transitions 

(qpts) i a a a a a km il u u m. 

Quantum phase transitions [lij are transitions between 
qualitatively distinct phases of quantum many-body sys- 
tems, driven by quantum fluctuations. In view of the 
connection between entanglement and quantum correla- 
tions 01 j one anticipates that entanglement will furnish 
a dramatic signature of the quantum critical point. Peo- 
ple hope that, by employing quantum entanglement, the 
global picture of the quantum many-body systems could 
be diagnosed, and one may obtain fresh insight into the 
quantum many-body problem. Hence, in addition to 
its intrinsic relevance with quantum information appli- 
cations, entanglement may also play an interesting role 
in the context of statistical mechanics. 

The aforementioned studies are based on the analysis 
of particular many-body models. Recently a theorem 
of the relation between QPTs and bipartite entangle- 
ment is proposed |l8j . The authors conclude that, under 
certain conditions, a discontinuity in or a divergence of 
the ground state concurrence [the first derivative of the 
ground state concurrence] is both necessary and sufficient 
to signal a first-order QPT (1QPT) [second-order QPT 
(2QPT)]. Most of the previous investigations for specific 
models support their conclusion. This may strengthen 
the belief that one can determine QPTs by using quan- 
tum entanglement. 

In this paper, the entanglement properties (the ground 
state concurrence and the von Neumann entropy) are cal- 
culated for an exactly solvable quantum spin model . 
Contrary to conventional wisdom, we find that there ex- 



ists a discontinuity in the first derivative of the concur- 
rence, at which there is no quantum critical point. In fact, 
similar result had already been discovered in Ref. 0] for 
a quantum spin model on a simplex in a magnetic field. 
Here we give general arguments to show why the an- 
alyticity property of the concurrence is more intricate 
than that of the ground state energy. Thus there is no 
one-to-one correspondence between QPTs and the non- 
analyticity property of the concurrence. That is, it is not 
always possible to infer the existence of QPTs from con- 
currence. Furthermore, for the one-dimensional XXZ 
model at the critical point of the isotropic ferromagnetic 
case, it is found that the first derivative of the concur- 
rence is discontinuous 0. However, it is a 1QPT, in- 
stead of 2QPT. The reason why the non-analyticity of 
the ground state energy of the XXZ spin chain does not 
faithfully reflect that of concurrence is explained. From 
Refs. [a,l2l an d our result, it is clear that QPTs in general 
can not be distinctly characterized through the analysis 
of the analyticity properties of concurrence. Moreover, 
we show that, for the model considered in this paper, 
the von Neumann entropy remains constant even cross- 
ing the critical point. That is, the von Neumann entropy 
can not always detect QPTs. Therefore, in order to have 
close connection with QPTs, some other measures of en- 
tanglement arc needed. 

The exactly solvable quantum spin model considered 
here is the isotropic spin-i XY (or spin-i XX ) chain 
with three-spin interactions [l9|. 

i=l 

+^K-i^JVi-^_i^i)], (i) 

where N is the number of sites, o~f (a = x, y, z) are 
the Pauli matrices, and A is a dimensionless parame- 
ter characterizing the three-spin interaction strength (in 
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unit of the nearest-neighbor exchange coupling) . Periodic 
boundary condition <tn+i = &i is assumed. This model 
can be solved by using the Jordan- Wigner transforma- 
tion [20I l2l| , and all physical quantities can in principle 
be calculated exactly. It is shown that the three-spin 
interaction can lead to a 2QPT at A c = 1 [l9| . 

Here we consider the entanglement between two spe- 
cific spins in the ground state of a quantum system. The 
state of two spins i and j in the ground state of a quantum 
system is described in terms of the two-particle reduced 
density matrix p^ obtained by tracing over other spins. 

The structure of pij follows from the symmetry prop- 
erties of the Hamiltonian. The Hamiltonian in Eq. 
is real and it has the following two symmetries. One is 
a global {/(l)-rotation symmetry about the spin- z axis, 
another is a Z2 symmetry of a global 7r-rotation about 
the spin-cc axis 12211. These symmetries guarantee that 
Pij has the form [23j 



with 
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in the standard basis {| TT>, I II), I IT), I II)}- Wang 
and Zanardi [24| have shown that the matrix elements of 
Pij can be expressed in terms of the various correlation 
functions (of o^) (a, = x, y, z): 



({of' 



(3) 



Note that 



> because of the inequality 



|(ofo|)| < 1, which is a special case of the Schwarz in- 
equality \{A*B)\ < vP^y^tB) for A = I (I is the 
identity operator) and B = afa^. 




uantify- 



From p^, the ground state concurrence |2 
ing the entanglement is readily obtained as |2 



Cij = 2 max {0, | Zy | - tiy } 

= \ max {0, 1 (of of ) + (o^)| - (a*a*) - l} (4) 

Because the entanglement between a pair of adjacent 
spins is expected to be dominant compared with a pair 
of non- nearest-neighbor spins, we focus on the nearest- 
neighbor case in the following discussions. 

Using the method adopted by Lieb, Schultz and Mat- 
tis |26| , one can calculate the spin-spin correlation func- 
tions exactly |Tgj. 
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We find that, for the nearest-neighbor cases, two correla- 
tion functions (of of +1 ) and (of of , x ) are dependent, and 
the latter can be written in terms of the former. Thus 
the nearest-neighbor concurrence is only determined by 
a single correlation function G. By substituting the re- 
sults of the correlation functions into Eq. Q , the exact 
expression of the concurrence between a pair of adjacent 
spins becomes 

C M+ i = imax{0,(G+l) 2 -2}. (7) 

The dependence of C^+i on A is plotted in Fig. 1. We 
see that, both at A = A c = 1 and A = Ao = 2/(v / 2— 1)tt ~ 
1.5369, the first derivative of the concurrence dC^i+x/dX 
shows discontinuities, while C^j+i is continuous. The 
discontinuity in dCi^+i/ d\ at A = A c = 1 do indicate 
the 2QPT of the present model, consistent with the pro- 
posal in Ref. 0. However, an unexpected discontinuity 
in dCi^+i/dX occurs at A = Ao, which is not a critical 
point! The origin of non-analyticity in the concurrence at 
A = Ao comes from the requirement that the concurrence 
should be non-negative, but not from the non-analyticity 
of p^. Therefore, the discontinuity in dCi^ + i/d\ needs 
not always indicate the existence of non-analyticity in the 
ground state energy and show any QPT. We note that 
the possibility of such an unanticipated discontinuity in 
dCi t i + i/dX is not addressed by the theorem of Ref. [lSj. 
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FIG. 1: The ground state concurrence of the nearest-neighbor 
spins d,i+i as a function of A for the XX chain with three- 
spin interactions in Eq. fl) . 



In a recent work studying the one-dimensional ex- 
tended Hubbard model |l5j, the authors show that QPTs 
can be identified at places where the von Neumann en- 
tropy is extremum or its derivative is singular. The von 
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Neumann entropy, another measure of entanglement, is 
defined as S = — tr(/jj log 2 pj), where pj is the one- 
particle reduced density matrix obtained by tracing over 
all sites except the j-th site, and therefore pj = trj(py ) 
(tr, stands for tracing over the i-th site). One may won- 
der if this measure of entanglement will still work for the 
present model. By using Eqs. J2J and 10 and after trac- 
ing over the i-th site for pij, one obtains pj = and 
the von Neumann entropy S = 1 for all A. The von Neu- 
mann entropy thus fails to detect the QPT of the present 
model. It is because the non-analyticity in the matrix el- 
ements of p^ cancel each other by taking trace over the 
i-th site. Hence the von Neumann entropy is not always 
useful to detect QPT. 

We note that the present example is not the only 
exception for the anticipation that the non-analyticity 
property of concurrence can be used to determine QPT. 
In general, the concurrence [2j| is defined by = 
max{0, Cij} with CV, = 71 - 72 - 73 - 74- Here j a 
are the four eigenvalues , in descending order, of the 



matrix Ri 



w = yvy y u i a '!i'>., a * ® a j)VPv Wltn Pij 
being the complex conjugate of the reduced density ma- 
trix p^ in the standard basis. We emphasize that, even 
for the models such that Cij is always non-negative and 
the matrix elements of p%j change smoothly as the phys- 
ical parameter, say A, is varied, the concurrence (in this 
case Cij — Cij) can still show a cusp-like singularity at 
some A = Ao- For example, this will happen (cf. Eq. (7) 
of Ref. ^ a t least one of the 7 Q takes the form of 

\A — B\ (where A and B denote two functions of A) and 
A — B changes sign at A = Ao- Therefore, the concurrence 
can be non-analytic, but it again does not correspond to 
a QPT. That is, although the eigenvalues j a are alge- 
braically related to the matrix elements of pij, Cij may 
still have different analyticity properties from p^. Re- 
minding that max{0, Cy } = \Cij\/2 + CV//2, it is clear 
that the two possibilities of the unexpected discontinu- 
ities in dC i>i+ i/dX discussed above all originate from the 
non-analyticity of the absolute-value function. As men- 
tioned before, an unexpected discontinuity in dCi^+i/dX, 
which does not indicate a QPT, had already been discov- 
ered in Ref. Q. We believe that the non-analyticity in 
the concurrence in that case may be due to the reason 
explained above. 

Even though the discontinuity in the first derivative 
of the concurrence does indicate a QPT, it may not be 
2QPT. An example is the one-dimensional XXZ model, 
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It is shown that, at the critical point A = — 1 (cor- 
responding to the ferromagnetic Hciscnberg model), 
dCi,i+i/9A is discontinuous, while C^+i is continu- 
ous and Ci i i + i\A=-i =0 0- However, it is a 1QPT 



(see below), instead of 2QPT. The reason why the non- 
analyticity of the concurrence of the XXZ spin chain 
does not faithfully correspond to that of the ground state 
energy is explained below. 

We first show the relations among the ground state 
energy (and its first derivative), the matrix elements of 
the reduced density matrix, and the concurrence for the 
XXZ spin chain. The concurrence of the XXZ spin 
chain has the same expression as Eq. (@J [3, lUJ, |2J, |2jj . 
Due to the translational invariance, the ground state 
energy per site for the XXZ spin chain can be writ- 
ten as £ = (afaf +1 ) + (cx?^) + M^fL-i)- Em- 
ploying the Hellmann-Feynman theorem |2g|. one has 
d£/dA = (of of +1 ). Thus, for the nearest-neighbor spins, 
two matrix elements Uj^+i and Zi i+i of the reduced den- 
sity matrix can be written as 



1 / d£ ' 
*.m 4 ( dA) ' 



(9) 



and the nearest-neighbor concurrence becomes C^+i 
max{0, Ci,i+i} with 



C, 



(£ + !) + (!- A) 



d£_ 
dA 



(10) 



Here we use the fact that {cfaf +1 ) + {o-\a^ +1 ) < 0. 
This inequality is satisfied for the XXZ spin chain in 
Eq. (JSJ, because the ground state wavefunction obeys 
the Marshall-Peierls sign rule [2^| . 

The XXZ spin chain is an exactly solvable model, and 
the expression of £ can be found in Ref. [3(J ■ For the crit- 
ical point A = — 1, one has £|a=-i = — 1> d£/dA — > 
as A -> -1+ and d£/8A = 1 for A < -1. Thus 
d£/dA is discontinuous at A = —1, which is a man- 
ifestation of a 1QPT. Based on these results, we find 
that C^j+i is indeed not continuous at A = — 1, where 
it has a finite jump from —1 to 0. However, because 
C^i+i = max{0, Cj^+i}, C^+i becomes continuous and 
equal to zero at A = — 1. That is, the discontinuity in 
Ci_i+i is hidden under the operation max{0, . . .}. That 
is the reason why the non-analyticity of the concurrence 
is not faithfully induced by that of the ground state en- 
ergy. In short, the discontinuity in the first derivative 
of £ (and therefore the matrix elements of the reduced 
density matrix) may not always lead to discontinuity in 
C^i+i. Therefore, a 1QPT may be misunderstood as a 
2QPT through analyzing the non-analyticity property of 
concurrence. 

There is another critical point of the XXZ spin chain 
at A = 1 (corresponding to the anti-ferromagnetic 
Heisenberg model). It is found that Cj^+i and 
dCi.i+i/dA are both continuous at A = 1, and 
<9CV;+i/(9A|a=i = (or Cj^+i reaches its maximum 
value at A = 1) 0, 0| . It is interesting to see how 
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these results can be realized in the present framework. 
At the critical point A = 1, it is shown in Ref. [3(J that 
£ and all of its derivatives with respect to A are contin- 
uous. Therefore, C^j+i (= C^j+i, because Cj.j+i > in 
this case) and SC^+i/dA will not show discontinuity at 
A = 1. Moreover, since 



t+i 



1 



dA 



(1-A) 



d 2 £ 
<9A 2 ' 



(11) 



we find that dCi.i+i/dA — > as A — > 1. Thus the results 
in Refs. 0,0 are reproduced. 

In summary, although many examples indicate that 
QPTs can be distinctly characterized through the ana- 
lyticity properties of concurrence, we stress in this pa- 
per that this viewpoint is not true in general. Except 
those cases of 2QPTs indicated by the discontinuity in 
<9C;,;+i/(9A, it is also possible that (i) dCij+i/dA is dis- 
continuous, but there is no QPT (Ref. and the present 
case); (ii) dCi^x/dA is discontinuous, while it is a 1QPT 
rather than 2QPT {XX Z spin chain at A = -1 0). 
We further point out that QPTs can not always be di- 
agnosed even by using the von Neumann entropy. As 
far as we know, there are some other measures of en- 
tanglement, say localizable entanglement [lfj and global 
measure of entanglement Therefore, while the ana- 
lyticity properties of concurrence and von Neumann en- 
tropy are not necessarily related to the existence of crit- 
ical points, other measures of entanglement may work. 
Thus more effort is necessary to clarify the relationship 
between QPTs and entanglement. 
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